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ABSTRACT 

Context. Hot temperature minerals have been detected in a large number of comets and were also identified in the samples of Comet 
Wild 2 that were returned by the Stardust mission. Meanwhile, observations of the distribution of hot minerals in yotmg stellar 
systems suggest that these materials were produced in the inner part of the primordial nebula and have been transported outward in 
the formation zone of comets. 

Aims. We investigate the possibility that photophoresis provides a viable mechanism to transport high-temperature materials from the 
inner solar system to the regions in which the comets were forming. 

Methods. We use a grid of time-dependent disk models of the solar nebula to quantify the distance range at which hot minerals can 
be transported from the inner part of the disk toward its outer regions as a function of their size (10"^ to IQ-' m) and density (500 and 
1000 kg m"^). These models will also yield information on the disk properties (radius of the inner gap, initial mass, and lifetime of the 
disk). The particles considered here are in the form of aggregates that presumably were assembled from hot mineral individual grains 
ranging down to submicron sizes and formed by condensation within the hottest portion of the solar nebula. Our particle-transport 
model includes the photophoresis, radiation pressure, and gas drag. 

Results. Depending on the postulated disk parameters and the density of particles, 10"^ to 10"' m aggregates can reach hehocentric 
distances up to ~35 AU in the primordial nebula over very short timescales (no more than a few hundred thousand years). 10"' m 
particles follow the same trajectory as the larger ones but their maximum migration distance does not exceed -26 AU and is reached 
at later epochs in the disks. On the other hand, 10"^ to 10 m aggregates are continuously pushed outward during the evolution of 
the solar nebula. Depending on the adopted disk parameters, these particles can reach the outer edge of the nebula well before its 
dissipation. 

Conclusions. Our simulations suggest that irrespective of the employed solar nebula model, photophoresis is a mechanism that can 
explain the presence of hot temperature minerals in the formation region of comets. Comets probably had the time to trap the dust 
transported from the inner solar system either in their interior during accretion or in the form of shells surrounding their surface if 
they ended their growth before the particles reached their formation location. 

Key words, planetary systems - protoplanetary disks - comets: general - comets: individual (8 IP/Wild 2) - Kuiper Belt: general - 
Oort Cloud 



1. Introduction 

Hot-temperature minerals have been detected in a large number 
of comets (Campins & Ryan 1989; Crovisier et al. 2000; Sitko 
et al. 2004; Wooden et al. 2000, 2004, 2010) and were also iden- 
tified in the samples of Comet 8 IP/Wild 2 that were returned 
by the Stardust mission (Brownlee et al. 2006). These miner- 
als include crystalline silicates that presumably condensed in 
the 1200-1400 K temperature range in the solar nebula (Hanner 
1999) and calcium, aluminum-rich inclusions (CAIs), which are 
the record of a very hot (1400-1500 K) stage of nebular evo- 
lution because they are composed of the first minerals to con- 
dense from a gas of solar composition (Grossman 1972; Jones 
et al. 2000). On the other hand, observations of young stellar 
systems show that the abundance of crystalUne silicates is much 
higher in the inner disk than in the outer disk, but that even the 
outer disks show more crystalline silicates than the interstellar 
medium (Tielens et al. 2005). These observations then suggest 



that crystalline silicates, and probably also CAIs, were produced 
in the iimer part of the primordial nebula and have been trans- 
ported outward in the formation zone of comets. 

A number of mechanisms has been invoked to account for 
the origin of these high-temperature minerals in comets. It has 
been proposed that shock waves in the outer solar nebula could 
anneal the amorphous silicates to crystallinity in situ prior to 
their incorporation in comets (Marker & Desch 2002). However, 
the isotopic composition, minor element composition, and even 
the range of Fe/Si ratios measured in the dust that was returned 
by the Stardust spacecraft from Comet 8 IP/Wild 2 appear to 
be inconsistent with an origin by annealing of interstellar sil- 
icates in the primordial nebula (Brownlee et al. 2006). An al- 
ternative possiblity is the radial mixing induced by turbulence 
which is responsible for the angular momentum transport within 
the primitive nebula (Shakura & Sunyaev 1973). This turbulence 
favors the rapid diffusion of the different gaseous compounds 
and gas-coupled solids throughout the nebula. One-dimensional 
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(vertically averaged) diffusive transport of particles in the disk 
(Bockelee-Morvan et al. 2002) or through its surrounding lay- 
ers (Ciesla 2007, 2009) has therefore been proposed to account 
for the presence of hot temperature minerals in the formation 
zone of comets. It is uncertain however whether turbulent trans- 
port suffices to explain the observations, or whether alternative 
physical processes are also needed. On the other hand, Hughes 
& Armitage (2010) recently studied the outward transport of 
particles in the nebula via a combination of advection (inward 
drift of particles though interaction with gas) and turbulent dif- 
fusion in an evolving disk. These authors found that the advec- 
tion of solids within the gas flow significantly reduces the out- 
ward transport efficiency for larger particles (typically a few mil- 
limeters), thereby limiting the extent of mixing uniformity that 
is achievable within the disk via turbulent diffusion. 

An alternative transport mechanism to turbulent diff'usion 
whose effects have been investigated in the last years in the solar 
nebula is photophoresis (Krauss & Wurm 2005; Wurm & Krauss 
2006; Krauss et al. 2007; Mousis et al. 2007; Wurm et al. 2010). 
This effect is based on a radiation-induced temperature gradient 
on the surface of a particle and the consequential nonuniform 
interaction with surrounding gas. When the existence of an in- 
ner gap is postulated in the disk, this latter becomes optically 
thin enough for particles to see the proto-Sun, but still has a rea- 
sonable gas content, which enables the photophoretic force to 
push dust grains outward (Mousis et al. 2007). This process pro- 
vides a mechanism to transport high-temperature material from 
the inner solar system to the regions in which the comets were 
forming. Eventually, the dust driven outward in this manner will 
reach a region where the gas pressure and irradiation are so low 
that the combined outward forces of radiation pressure and pho- 
tophoresis can only balance the inward drift of particles. 

In this work, we use a grid of time dependent models of the 
solar nebula to quantify the distance range at which particles (i.e 
hot minerals) can be transported from the inner part of the disk 
toward its outer regions as a function of their size and density 
as well as of the disk properties (radius of the inner gap, initial 
mass, and lifetime of the disk). The grid of models used here 
allows us to consider the full range of thermodynamic condi- 
tions that might have taken place during the solar nebula's evo- 
lution. The particles considered in our model are in the form 
of hot mineral aggregates with sizes ranging between 10"^ and 
10"' m. The trajectories of particles with lower sizes are gen- 
erally influenced by radiation pressure while those of particles 
with larger sizes begin to be mostly affected by gas drag. The ag- 
gregates are presumed to have been assembled from hot mineral 
individual grains ranging down to submicron sizes. We consider 
that these hot minerals have formed by condensation within the 
hottest portion of the solar nebula, well inside 1 AU (Chick & 
Cassen 1997). We also show that the determination of the dust 
size distribution within rings observed in young circumstellar 
disks and their position relative to the parent star is likely to 
bring some constraint on the lifetime and eventually the initial 
mass of the disk from which they originate. 

Section 2 is devoted to the description of our modeling ap- 
proach, detailing the particle transport and solar nebula models 
employed in this work. In Section 3 we detail the disk and par- 
ticle parameters employed in our different models. In Section 4 
we present and analyze the trajectories of particles determined in 
the frame of these models. In Section 5 we show that calculations 
of particle trajectories induced by photophoresis can be used as 
a tool to determine some physical parameters of circumstellar 
disks. Section 6 is devoted to the discussion of the assumptions 
of our model. 



2. Model 

2.1. The photophoretic force 

Any particle embedded in gas and heterogeneously heated by 
light feels a photophoretic force, which usually pushes it away 
from the light source (Krauss & Wurm 2005; Wurm & Krauss 
2006). The force is strongly pressure-dependent and can be 
stronger than radiation pressure and the gravity of the Sun by 
orders of magnitude in the solar nebula. This mechanism in- 
duces the migration of particles ranging from micron to cen- 
timeter sized in the solar nebula under the combined action of 
photophoresis, radiation pressure, and gas drag, provided that 
the disk is sufficiently transparent (Mousis et al. 2007). 

Following the approach developed by Krauss et al. (2007), 
we assume here that the disk's gas flow conditions are described 
by the Knudsen number, Kn, which is defined as Kn - 1 1 a, 
where I is the mean free path of the gas molecules and a is the ra- 
dius of the particle. If the mean free path of the gas molecules is 
large compared to the considered particle sizes, i.e., for Kn > 1, 
then the gas flow is in the free molecular flow regime. In the con- 
trary case, the gas flow is in the continuum regime. In these con- 
ditions, the photophoretic force Fp/, on a spherical particle, valid 
for both flow regimes, can be expressed as follows (Beresnev et 
al. 1993): 



Fph = -fl /7i I 



\2kTl aE + l5AKn{\-aE)l^ + aEK^2' 



(1) 



where / is the light flux (power incident per area), nig is the av- 
erage mass of the gas molecules (3.89 x lO^^^kg), T is the gas 
temperature, and k the Boltzmann constant. J\ is the asymmetry 
factor that contains the relevant information on the distribution 
of heat sources over the particle's surface upon irradiation. In the 
following calculations, we assume J\ - 0.5, which corresponds 
to the case where the incident light is completely absorbed on the 
illuminated side of the particle. The energy accommodation co- 
efficient ue is the fraction of incident gas molecules that accom- 
modate to the local temperature on the particle surface and, thus, 
contribute to the photophoretic effect. Here, we assume complete 
accommodation, i.e, ue - 1- 

The thermal relaxation properties of the particle are sum- 
marized in the heat exchange parameter A = A^ff/Ag, where 
Ag is the thermal conductivity of the gas and A^ff the effective 
thermal conductivity of the particle. For Ag, we adopt values for 
molecular hydrogen for temperatures above 150K (as tabulated 
by Incropera & DeWitt 2002). On the other hand, because he- 
lium has a higher thermal conductivity than hydrogen for lower 
temperatures, we assume that this species determines the ther- 
mal conductivity of the gas and use values taken from the com- 
pilation of Bich et al. (1990). The expression of A includes the 
conduction of heat through the particle and the thermal emission 
from the particle's surface, according to 



A^ff - Ap + AecrT a. 



(2) 



where Ap is the thermal conductivity of the particle supposed 
here to be 10"^ Wm"' K"' (Mousis et al. 2007), e its emissivity 
assumed to be 1, and cr the Stefan-Boltzmann constant. On the 
other hand, the functions ijji and 1/^2 in Eq. ([TJ depend only on 
Kn in the form 



Kn 

Kn+(5nl\i) 



(1 



^2 = (i + f^«)(l 



1.2In 



-Kn 



W0Kn-+7r/4 



)■ 



(3) 
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As noted by Krauss et al. (2007), an additional photophoretic 
force arises if the accommodation coefficients vary over the sur- 
face of the dust grain (Cheremisin et al. 2005), but we restrict 
the treatment to the "classical" photophoretic force as given in 
Eq. ([T]). In the present work, we assume that Eq. ([T]i is valid for 
all parts of the solar nebula and all particles. 



2.2. Balistic transport 

In a protoplanetary disk where the gas pressure (in the midplane) 
decreases with distance from the star, the gas is supported by a 
pressure gradient and rotates slower than the Keplerian veloc- 
ity (Weidenschilling 1977). Solid particles are only stable on a 
Keplerian orbit. Therefore, interaction with the gas leads to an 
inward drift of solids toward the star For particles that couple 
to the gas flow on timescales short compared to an orbital pe- 
riod, the problem reduces to a one-dimensional (radial) calcula- 
tion. The inward drift is then induced by the fraction of gravity 
(residual gravity), which is not balanced by the circular motion 
with the sub-Keplerian gas velocity. The force, Fies, acting on a 
particle of mass nip due to residual gravity is given as 



nip dp 
Pg dr' 



(4) 



where pg is the density and p the pressure of the gas. 

In addition, radiation pressure has also to be considered for 
at least micron-sized particles (Krauss & Wurm 2005). The ra- 
diation pressure force can then be expressed as follow; 



Frad = TTfl 



C light 



(5) 



where cugiu is the speed of light. The sum of the outward forces 
(Eq. ([T]) and Eq. (j5]l) and the inward force (Eq. (j4|) gives the drift 
force Fdrifr- We treat the problem as being purely radial here be- 
cause we are mostly interested in the small particles. These small 
particles couple to the gas on timescales much shorter than the 
orbital timescale, which justifies the radial treatment as outlined 
in Wurm & Krauss (2006). The radial drift velocity with respect 
to the nebula is then estimated to be 



^ ph ^ md "t" ^ res 
Vdr = T, 



(6) 



where r is the gas grain coupling time and nip the mass of the 
considered particle. 

As larger dust aggregates drift outward, they pass from a re- 
gion where the continuum flow regime is valid to a region where 
the free molecular flow regime applies. Hence, as with the pho- 
tophoretic force, we have to consider an equation describing the 
gas grain friction time in both regimes. It is given by 



6nria 



C,, 



(7) 



where 77 is the dynamic viscosity of the gas. This assumes Stokes 
friction, which is justified because the Reynolds numbers for 
the drift of particles smaller than 10 cm are well below 1. The 
Cunningham correction factor, Cc, accounts for the transition 
between the different flow regimes (Cunningham 1910) and is 
given as (Hutchins 1995) 



Cc = 1 (1.231 +0.47e 



-1.178/A:n\ 



(8) 



To close the set of equations, we need to determine the dy- 
namic viscosity rj and the mean free path I. In the framework of 
the classical kinetic theory for dilute gases (see e.g. Reif 1972), 
these quantities are given by 



1 



8kT . 



and 



^/2n 



(9) 



(10) 



where n is the molecule number density in the gas, and cr the 
collisional cross section of the gas molecules. The latter is very 
difficult to obtain. It is easier to find the value of the dynamic vis- 
cosity at a given temperature for H2 and then use the functional 
form of 77 in Eq. |9]to determine its value at any temperature. We 
use 770 = 9.0 X 10 Pas at To = 300 K (Lide 2007). Finally, 
inverting Eq. |9] one obtains the mean free path. 

We note in passing that Beresnev et al. (1993) used a nor- 
malising factor of 1/2 in Eq.|9]instead of the 1/3 that applies for 
three-dimensional gases, and this may slightly modify the nu- 
merical constants in Eq.[T] However, the change is likely smaller 
than the uncertainties because of all the approximations made to 
solve the conservation equations in their model. 

Our description of the radial transport of particles in the disk 
includes their drag back toward the central star by the infalling 
nebula flow that moves at the velocity of Vac- In our disk model, 
the accretion speed ranges from a few tens of cm/s in the inner 
part to below one cm/s at larger distance in the early stages and 
substantially decreases later on. In the simplified solar nebula 
model presented in Section 2.3 the accretion velocity is esti- 
mated to be 



2t ■ 



(11) 



where r is the distance from the Sun, 



tvis = 3^ ^ ^ is the typical 
local viscous time, H is the local height of the nebula, Q is the 
local Keplerian frequency given by = GMqIP and a is the 
viscosity parameter of the disk described in Section [23] Finally, 
the position of particles is integrated from the inner edge of the 
disk at time f = to a position r{tdisk) at the age of the disk t^isk 
via 



r(tdisk) 



-I ' 



(VdMt),t)-VaMt),t))dt. 



(12) 



2.3. The protoplanetary disk 



The structure and evolution of the protoplanetary disk is mod- 
eled as a non irradiated, 1 + lD turbulent disk, following the 
method originally presented in Papaloizou & Terquem (1999) 
and also developed by Alibert et al. (2005). The diffusion equa- 
tion (see Lynden-Befl & Pringle 1974; Papaloizou & Lin 1995) 
describing the evolution of the gas surface density E is conse- 
quently solved as a function of time t and distance r to the star; 



dl. 
dt 



3d_ 
r dr 



dr 



(13) 



where E is the surface density of mass in the gas phase in the 
nebula and v the mean (vertically averaged) turbulent viscos- 
ity. Compared to the original equation, the photo-evaporation 
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term, £,v(r), was added and is taken to be the same as in Veras 
& Armitage (2004). The mean turbulent viscosity is determined 
from the calculation of the vertical structure of the nebula: for 
each radius, r, the vertical structure is calculated by solving 
the equation for hydrostatic equilibrium together with the en- 
ergy equation and the diffusion equation for the radiative flux 
(see Papaloizou & Terquem 1999). The local turbulent viscosity 
(as opposed to that averaged in the vertical direction) is com- 
puted using the standard Shakura & Sunyaev (1973) formalism: 
V = aCllQ., where is a free parameter and C, the local speed 
of sound determined by the equation of state. Using this proce- 
dure, we derived the midplane pressure and temperature as well 
as the mean turbulent viscosity as a function of r and S. These 
laws are finally used to solve the diffusion equation (Eq. (13 i) 
and to calculate the pressure- and temperature-dependant forces 
on dust grains. Figure 1 represents the temperature, pressure and 
surface density profiles in the midplane of the disk characterized 
by a mass of 0.03 Mq and a lifetime of 6 Myr (see Sec. 3 for fore 
details) at different epochs of its evolution. 

Following the approach of Mousis et al. (2007), we con- 
sider that the disk is not optically thin and that Rayleigh scat- 
tering from molecular hydrogen is the dominant dimming effect 
in the nebula (Mayer & Duschl 2005) for temperatures below 
1500 K and at wavelengths shorter than a few jim. This con- 
dition is fulfilled only after 10^ yr and beyond 0.5 AU in all 
the solar nebula models used in our calculations. For H2, i.e, 
the dominant molecule, the Rayleigh scattering cross section 
is cr(X) = 8.49 X IQ-'^^ I A'^icrn^) (Vardya 1962). Assuming the 
illuminating light follows a black body spectrum, the Planck 
mean cross section as a function of the black body temper- 
ature Tb is found to be a-{TB) = 1.54 x 10""*^ (cm^) 
(Dalgarno & Williams 1962). Note that, in our case, Tb is not 
the temperature of the nebula, but rather the effective temper- 
ature of the illuminating source, the Sun. With a disk's mean 
molar mass of 2.34 g/mol, the mass absorption coefficient is 
found to be a-,„{TB) = 3.96 x 10"'^ (cm^/g). The effec- 
tive temperature and the luminosity of the early Sun were taken 
from the ZAMS (Zero Age Main Sequence) model computed 
by Pietrinfemi et al. (2004), which is available in the BaSTI 
database (http://albione.oa-teramo.inaf.it). We chose the param- 
eters relevant for the Sun, i.e, a solar mixture of heavy ele- 
ments, no overshooting, a metallicity Z- 0.0198, and a helium 
content Y- 0.273 (Z and Y together in the mass fraction). In 
this model, the surface temperature of the early Sun is 5652 
K and its initial luminosity is 2.716 x 10^'' W. We derived 
cr„(5770) - 4.0 X lO"'* (cm^/g) from the adopted effective tem- 
perature of the early Sun. Light becomes extinguished close to 
the star as a result of the high gas density, while the outer regions 
play only a minor role in the extinction. 

3. Choice of parameters 

We constructed a grid of nine disk models encompassing the 
range of thermodynamic conditions that might have taken place 
during the solar nebula's evolution. The three initial disk masses 
were fixed to 0.01, 0.03 and 0.1 Mq respectively, with 0.01 
Mq corresponding to the minimum mass solar nebula (hereafter 
MMSN) defined by Hayashi (1981). The initial mass of each 
disk is integrated between 0.25 and 50 AU and the initial gas 
surface density is given by a power law Z oc r^^l^, with an ini- 
tial value taken to be Z(5.2AU) = 100, 300, and 1000 gcm^^ at 
5.2 AU for disk masses of 0.01, 0.03 and 0.1 Mq, respectively. 
Here, the lifetime of the disk is governed both by viscosity and 
photoevaporation by the Sun or nearby stars. On the other hand. 



the viscosity parameter rules the accretion velocity of the disk 
(Eq. 12) but this latter is found to be low compared to the veloc- 
ities due to photophoresis and gas drag for particles larger than 
10""* m (see Fig. 2 for an example of particle velocities due to 
photophoresis, radiation pressure, residual gravity and accretion 
flow along their trajectories in the nebula). Here the viscosity 
parameter is fixed to 7 x lO""*, i.e, a value adopted in works aim- 
ing at synthesizing different populations of planets around other 
stars (Mordasini et al. 2009a, 2009b) and the photoevaporation 
rate is varied to obtain the appropriate disk lifetimes (1,3, and 
6 Myr for each selected mass). In each case, the lifetime coiTe- 
sponds to the time taken for the mass of the disk (integrated until 
50 AU) to decrease to 1 % of its initial value. 

Mousis et al. (2007) have calculated the optical depth of the 
disk at 30 AU as a function of time. They found that even at 
late epochs, only ~0.1% of the Sun's radiation is available in 
this region. As a result, these authors found that the high ex- 
tinction induced by H2 Rayleigh scattering limits the outward 
transport of particles only to very short heliocentric distances 
(typically a few AU) when they are released from the inner- 
most regions. On the other hand, particle transport can be en- 
hanced at larger heliocentric distances when a gap is formed in 
the inner disk. In particular, there is a growing body of obser- 
vational evidence for the existence of disks whose inner few AU 
are cleared or are strongly depleted of gas (D' Alessio et al. 2005; 
Sicilia-Aguilar et al. 2006; Espaillat et al. 2008; Pontoppidan et 
al. 2008; Thalmann et al. 2010). For this reason, and similar to 
Mousis et al. (2007), we assume here the presence of 1 and 2 
AU inner gaps within the nebula during the course of its viscous 
evolution. Gaps are prescribed in a way independent of the struc- 
ture of the disk models used in this work and their sizes remain 
constant with time. As shown in Sect. |4) such an inner hole is 
large enough to leave a reasonable fraction of the incoming light 
to let photophoresis work even in the outer solar system. 

Particles considered in our simulations have sizes ranging 
between 10"^ and 10^'m and are assumed to be spherical and 
composed of olivine, with a variable porosity. Density of aggre- 
gates is varied between 500 and 1000 kg m"^. The first value cor- 
responds to the random deposition of irregular olivine particles 
with density of 3300 kgm^^*, with a 15% filling factor (Blum 
& Schrapler 2004). The second value corresponds to the aver- 
age density measured in cometary interplanetary dust particles 
(Joswiak et al. 2007). We do not consider particles with sizes 
lower than 10"^ m because their path in the nebula is essentially 
controlled by radiation pressure. Moreover, for objects larger 
than about 1 m, the radial treatment we apply does no longer 
hold because the gas grain friction times become comparable to 
the orbital period. 

4. Outward transport of hot temperature aggregates 

All our calculations are based on the assumption that the disk 
opacity is essentially caused by Rayleigh scattering and not to 
dust, implying that the dust size distribution in the nebula is 
dominated by large particles instead of small particles. In the 
contrary case, smallest aggregates (here 10"^ m) would create 
a prominent opacity in the disk, implying that larger aggregates 
could only migrate outward in the wake of the small ones. 

Figures 3-6 represent the trajectories of 10"^ to 10"' m ag- 
gregates in the solar nebula that were computed using the defined 
particle densities and a set of six disk models that are expected 
to encompass the range of plausible thermodynamic conditions 
within the solar nebula (disk masses of 1 MMSN, 3 MMSN, 10 
MMSN with lifetimes of 1 or 6 Myr). At the beginning of each 
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Fig. 2. Velocities of particles due to photophoresis (a), radiation pressure (b), residual gravity (c), and accretion flow (d) represented 
as a function of time and size in the midplane of the disk characterized by a mass of 0.03 and a lifetime of 6 Myr. The density 
of particles is 500 kg m"^ and the radius of the inner gap is 2 AU. Position of larger particles essentially corresponds to the balance 
between photophoresis and residual gravity velocities. When the disk opacity is prominent, i.e., at early epochs, the position of 
smaller particles is mainly driven by the balance between photophoresis and accretion flow velocities. At later epochs, the position 
of these particles becomes ruled by the balance between all velocities. 



computation, the particles start their migration within the disk 
from the outer edge of the inner gap. 

Figure 3 shows that 10^^-10"' m particles with densities of 
500 kg m"^ that migrate within a disk with a 1 AU inner gap 
can reach heliocentric distances ranging between ~24 and 28.1 
AU, depending on the choice of the initial mass and lifetime of 
the nebula. Each of these positions corresponds to an equilibrium 
reached at the position where the outward drift of aggregates just 
balances the accretion flow and in no more than a few hundred 
thousand years. With time, the location of these particles slightly 
rebounds toward the Sun until the dissipation of the disk. The 
figure also shows that 10"^ m particles follow the same trajectory 
as the larger ones but their equilibrium position is reached at 
lower heliocentric distance (~20-22.8 AU) and at later epochs 
in disks owning similar input parameters. Interestingly enough, 
the position of smaller aggregates (10^^-10""* m) continuously 
progresses outward during the evolution of the disks. 10"^ m 
particles can even be pushed beyond the outer edge (~50 AU) of 



the nebula if one selects a low-mass disk (1 MMSN) with a long 
lifetime (6 Myr). This is because of the strong decrease of the 
gas density and opacity in this model that enables the radiation 
pressure to push the particles at higher heliocentric distance. 

Figure 4 represents the trajectories of the same particles 
as in Fig. 3, but for disk models with inner gaps fixed to 2 
AU. Because the Rayleigh scattering through H2 is strongly di- 
minished here, all particles reach higher heliocentric distances 
than in the cases considered in Fig. 3, but for similar migration 
timescales. Thus, 10^^-10 ' m particles reach heliocentric dis- 
tances as high as ~27-35.1 AU, depending on the adopted pa- 
rameters of the disk. In similar conditions, 10"^ m particles are 
also able to reach the ~23.9-26.3 AU distance range within the 
nebula. Moreover, 10"^ and even 10""^ m particles reach the edge 
of the nebula for low mass (1 MMSN) and long Ufetime (6 Myr) 
disk. 

Figures 5 and 6 show the trajectories of 10"^ to 10"' m ag- 
gregates with densities of 1000 kg m"^ within disk models with 
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inner gaps of 1 and 2 AU, respectively. Migration timescales 
remain similar to the previous cases: larger particles migrate 
very rapidly toward a maximum heliocentric distance while 
smaller ones continuously drift outward during the evolution 
of the disk. Because (i) the inward drift linearly depends on 
the mass of the aggregate (see Eq. 4) and (ii) the radial drift 
velocity v^r is inversely proportional to this quantity (see Eq. 6), 
all particles here migrate at lower heUocentric distances than in 
cases of disks based on similar parameters. Indeed, 10^^-10"' 
m particles do not exceed -20.2-23.6 AU (-22.7-29.7 AU) in 
the case of disks with 1 AU (2 AU) inner gaps. The maximum 
migration distance reached by intermediary size particles (10"^ 
m) becomes ~20 AU (~23.3 AU) in the case of disks with 1 
AU (2 AU) inner gaps. In every case, the maximum migration 
distance of lO""* m particles is several AU smaller than those of 
same size particles with densities of 500 kgm~^. Now only 10"^ 
m particles reach the outer edge of the nebula for low mass (1 
MMSN) and long lifetime (6 Myr) disk, irrespective of the gap 
size. 



5. Probing the dissipation of circumstellar disl(s 

Particle transport through the combination of photophoresis and 
radiation pressure has been invoked to explain the presence of 
ring-shaped dust distributions in young circumstellar disks such 
as the one around HR 4796A (Krauss & Wurm 2005). Here we 
show that in some cases, the determination of the dust size distri- 
bution within rings and their position relative to the parent star 
is likely to bring some constraints on the lifetime and eventu- 
ally on the initial mass of the circumplanetary disk from which 
they originate. Indeed, Figure 7 represents the settling distances 
reached by particles of different sizes and with densities of 500 
kg m"^ at the end of the solar nebula evolution. The figure shows, 
for example, that ring-like structures essentially composed of 
IQ-^ m particles and located in the -1.8-12 AU (3-21 AU) dis- 
tance range from the star could have formed in disks with 1 AU 
(2 AU) inner gap, which have short or intermediary lifetimes 
(here 1-3 Myr), irrespective of the initial disk's mass. In addi- 
tion, same size particles located at long distance to the star, i.e., 
-50 AU (upper hmit owing to the truncation of our model) or 
farther, could have formed in disks with long Ufetimes (6 Myr) 
and low initial masses (1 MMSN), irrespective of the size of the 
inner gap. To a lesser extent, one can also identify in Fig. 7 a 
relationship between the position (in the 15-25 AU range) of 
ring-hke structures dominated by the settUng of 10"^ to 10"^ m 
particles and the disk's Ufetime and irmer gap size. 

6. Discussion 

6.1. Disk's structure and evolution 

One could argue that the existence of an inner gap at early epochs 
within the nebula remains questionable. Indeed, gaps are often 
found in disks (i.e., transition disks) that are millions of years 
old. In this context, a significant offset might exist between the 
times of the difl'erent models used in this work and the chronol- 
ogy of the solar system formation and evolution that is testified 
by meteorite measurements or by the age of disks as inferred 
from luminosity studies of protostars. For these reasons, the less 
massive disk models used in this work that are associated to in- 
ner gaps correspond to cases that are the most consistent with 
the structure of transition disks. Moreover, our calculations are 
based on the assumption that the nebula is essentially devoid of 



dust, i.e., that the dust opacity is negligible. Indeed, if we assume 
that the smallest aggregates (10^^ m) have created a prominent 
opacity in the nebula, larger aggregates could only follow the 
small ones and reach the formation zone of comets toward the 
end of the disk evolution. 

Despite these caveats, the use of a set of disk models cov- 
ering the whole range of plausible thermodynamic conditions 
that took place in the primordial nebula allows us to show that 
hot-temperature minerals can drift up to heliocentric distances 
reaching -34 AU for the largest particles and 50 AU or beyond 
for the smallest ones provided that i) the existence of an in- 
ner gap is postulated within the nebula and ii) the opacity of 
the smallest dust particles remains neghgible inside the pho- 
tophoretic transport front. These simulations suggest that, irre- 
spective of the employed solar nebula model, photophoresis is 
a mechanism that can explain the presence of hot-temperature 
minerals at early epochs of the disk's evolution in the formation 
region of comets (from 10 to 30 AU according to the difi'erent 
scenarios - see, e.g.. Homer et al. 2007 for a review). Because 
comets have presumably accreted within a few hundred thou- 
sand years (Weidenschilling 1997), i.e., a timescale shorter than 
the one probably required to form an inner gap in the disk and to 
allow the photophoretic transport of particles formed close to the 
Sun, they probably had the time to essentially trap the dust trans- 
ported from the inner solar system in the form of shell surround- 
ing their surface if their accretion ended before particles reached 
their formation location. Because photophoresis works hetero- 
geneously, depending on the individual properties of a dust ag- 
gregate (composition, size, thermal and optical properties), one 
would expect the bulk of comets to be laden with particles of 
size and/or composition that would vary as a function of their 
accretion distance in the solar nebula. It is important to note that 
our calculations were made with the typical values of 0.5 and 
1 for the asymmetry factor Ji and the emissivity e. Assuming 
lower but still reasonable values for these parameters would not 
alter our conclusions. For example, if one assumes Ji = 0.4 and 
e = 0.8 in the disk model characterized by a mass of 3 MMSN 
and a lifetime of 6 Myr, 10 '-10^^ m particles still migrate up 
to - 25.4-25.9 AU (29.5-30 AU) and lO-^-lO""* m particles up 
to - 1 1-23.8 AU (18.5-33.3 AU) in the nebula owning a 1 AU 
(2 AU) inner gap. 

6.2. Role of turbulence 

The influence of turbulence on the particle motion has also to be 
considered in comparison to photophoresis. Two diff'erent cases 
have to be discussed. In the first case, the disk has an inner 
clear region and an optically thick outer region where opacity 
is provided by dust grains. In principle photophoresis is capable 
of moving the edge between optically thin and thick parts out- 
ward, thus clearing the disk inside-out from solids (see Krauss 
et al. 2007 for details). In this scenario, turbulent inward diffu- 
sion might counteract the outward motion of the edge. However, 
Krauss et al. (2007) discuss that turbulence will not prevent the 
outward motion of the edge. The edge wiU finally reach a lower 
heUocentric distance than in a similar situation where turbulence 
is negligible. The position of the edge will depend on the effec- 
tive inward transport of particles by the turbulence. The second 
case would be more Uke the situation discussed in this paper. In a 
disk where solid particles are treated as test particles and opacity 
is essentially generated by the gas (Rayleigh scattering) and not 
by dust extinction, turbulence is not an issue. It only broadens 
any ring-like particle concentration because photophoresis is a 
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Fig. 3. Position of particles of size 10 ^ to 10 ' m, as a function of time for disks with masses of 1, 3, or 10 MMSN and lifetimes of 
1 or 6 Myr. The density of particles is 500 kgm"^ and the radius of the irmer gap is 1 AU. 




Fig. 4. Same as in Fig. 3, but for an inner gap radius of 2 AU. 
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Fig. 5. Position of particles of size 10 ^ to 10 ' m, as a function of time for disks with masses of 1, 3, or 10 MMSN and lifetimes of 
1 or 6 Myr. The density of particles is 1000 kgm"^ and the radius of the irmer gap is 1 AU. 




Fig. 6. Same as in Fig. 5 but for an inner gap radius of 2 AU. 



Moudens et al.: Photophoretic transport of hot minerals in the solar nebula 



9 



directed force while turbulence is diffusive and acts statistically 
in both directions. 

6.3. Influence of particle rotation 

Particle rotation is an interesting topic for photophoretic forces. 

Photophoresis depends on the fact that a temperature gradient is 
established across the particle. Usually, an illuminated particle 
is warmer on the bright side than on the dark side. This temper- 
ature gradient always needs a certain time to adjust to changes 
in the illumination. Therefore, if the particle rotates, the temper- 
ature gradient might be different from the simple case where the 
particle is considered at rest. Rotations might be divided into two 
classes. 

The first class corresponds to rotations around the direction 
of the incident radiation. This rotation does not change the front 
and back side of the particle. The temperature gradient in a refer- 
ence frame fixed to the particle always stays the same. While the 
main component of photophoresis will still point away from the 
light source, the subtleties of particle morphology and composi- 
tion will induce sideward components of photophoresis as well. 
If a particle only rotates around the direction of light, the side- 
ward motion will oscillate but the component along the direction 
of Ught will remain constant. 

The second class corresponds to rotations around axes per- 
pendicular to the direction of light. These rotations change front 
and back of the particle. For slow rotations, the temperature gra- 
dient can follow the rotation. For somewhat faster spin the gra- 
dient lags behind, meaning that the cold side will trail into the 
warm side and vice versa. In total this results in a certain tem- 
perature gradient directed in a direction perpendicular to the di- 
rection of Ught. On a given orbit, these directed forces along 
the orbit can accelerate or slow down the particle and make it 
drif inward or outward. This is a photophoretic analog to the 
Yarkovsky effect, which works by radiation pressure, but the 
photophoretic effect can be several orders of magnitudes larger. 
This rotation would change the calculations given here one way 
or the other. Wurm et al. (2010) present first experimental evi- 
dence showing this eflfect in microgravity experiments for rotat- 
ing particles. 

However, this second class of rotations can only occur tem- 
porarily and is not of general importance. To see this, we remind 
the reader that the idea of photophoresis is founded on the fact 
that particles in protoplanetary disks are embedded in a gaseous 
environment. This point might explicitly be noted in view of par- 
ticle rotation. This simple fact is important because it imphes a 
fundamental difl'erence to particles in the current solar system 
where particles move in vacuum. While a particle in the solar 
system retains any arbitrary rotation state for a long time, ran- 
dom rotation of small particles in protoplanetary disks is rapidly 
damped away. For example, if a particle rotates through colh- 
sions with other particles, this rotation is gone on the order of 
one gas-grain friction time, which for bodies smaller than a me- 
ter is much shorter than the orbital timescale. Typical values for 
dust particles in the inner disk would be on the order of sec- 
onds, depending on the model and specific location in detail. 
Therefore, to retain a rotation, a constant torque around a given 
axis has to be applied to the particle. There are only two effects 
that can lead to torques: 

1. Particle motion with respect to the gas, i.e., radial drift. In 
analogy to particles in Earth's atmosphere (like snowflakes), 
particle rotation will be around the drift axis. For the small 
particles considered here, this relative motion can be re- 



garded as purely radial and transversal motion is not impor- 
tant (Weidenschilling 1977). A windmill might be an appro- 
priate visualization for this. 
2. Radiation-induced torques. In principle this might be re- 
garded in analogy to gas drag (assuming photons instead of 
gas molecules) and any rotation is around the direction of 
illumination. 

In both cases, a torque around another axis might initially be 
present but this will align the particles and only the systematic 
torques around the drift direction or radiation direction pertain. 
In both cases a potential rotation axis is oriented toward the star. 
This will always adjust during the orbit. This rotation does not 
change the front and back side and does not influence radial pho- 
tophoresis. Details with quantitative estimates can be found in 
Krauss et al. (2007). 

6.4. Prospects 

A more reahstic description of the disk's structure requires one 
to account for irradiation by the Sun at the disk atmosphere's 
surface (D'Alessio et al. 1998; Hueso & Guillot 2005; Garaud 
& Lin 2007; Cabral et al. 2010). Comparisons between irradi- 
ated and non irradiated models show that, for similar disk pa- 
rameters, the midplane temperature in the outer part of the disk 
becomes substantially higher (up to a few dozen of K) in the first 
series of cases (Garaud & Lin 2007; Cabral et al. 2010). In order 
to estimate the influence of irradiation on photophoretic trans- 
port, we need to know not only the temperature profile, but also 
the pressure and volume density in the midplane of the disk. As 
a first attempt, we increased the temperature by 30 K or up to 
100 K for any initial temperature lower than 100 K in our nom- 
inal disk model, to mimic the difference in temperature caused 
by irradiation. The trajectories of larger particles (10~^-10"^ m) 
remain almost similar to those plotted in Figs. 3-6, except for 
the maximum distances that are 2 to 5 AU closer to the Sun. The 
trajectories of small particles (10^^ m) are more affected by the 
temperature difference because their maximum migration dis- 
tance is about half that shown in Figs. 3-6. Modifying the mid- 
plane density to keep the product Txp constant roughly doubles 
the effect. However, both the photophoretic force and gas drag 
depend on the gas pressure and density and a fully consistent 
irradiated disk model will be needed to investigate the real influ- 
ence of irradiation on the trajectories of transported particles. 

An interesting evolution of this work would also be to con- 
sider explicitly the motion of dust particles in both radial and 
vertical directions. Indeed, the disk's upper layers are more 
transparent than those close to the midplane (because of dust 
sedimentation), and are therefore a perfect place for photophore- 
sis to be effective in the earliest phase of the disk's evolution 
and prior to the "transition disk" phase. At this epoch turbulence 
driven by active magnetorotational instability (MRI) may even 
be a useful ingredient because it would help a fraction of parti- 
cles to be maintained above the photosphere and be transported 
outward thanks to both photophoeresis and radiation pressure. 
Indeed, computations of turbulence in MRI disks (Turner et al. 
2010) have shown that turbulence is increasingly effective with 
scale-height, and consequently, may help to maintain particles 
high above. This may be a potentially interesting mechanism 
that would help photoporesis to be effective even in the youngest 
ages of the disk. 
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Fig. 7. Heliocentric distance reached by aggregates at the end of 
the solar nebula evolution as a function of the disk parameters 
and of the particle sizes. Density of all particles is 500 kg m 
and the inner gap of the disk is 1 AU (top) and 2 AU (bottom). 
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Fig. 1. Temperature, pressure and surface density profiles in the 
midplane of the disk characterized by a mass of 0.03 Mq and a 
lifetime of 6 Myr From top to bottom in each panel, times are 1, 
2, 3, 4, 5, and 6 Myr. 



